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' Abstract 

o : 

For a large class of field theories there exist portions of parameter 
space for which the loop expansion predicts increased symmetry break- 
ing at high temperature. Even though this behavior would clearly have 
far reaching implications for cosmology such theories have not been 
fully investigated in the literature. This is at least partially due to the 
pi il counter intuitive nature of the result, which has led to speculations 

^ ' that it is merely an artifact of perturbation theory. To address this 

issue we study the simplest model displaying high temperature sym- 
metry breaking using a Wilson renormalization group approach. We 
^ ■ find that although the critical temperature is not reliably estimated by 

the loop expansion the total volume of parameter space which leads 
to the inverse phase structure is not significantly different from the 
perturbative prediction. We also investigate the temperature depen- 
dence of the coupling constants and find that they run approximately 
according to their one-loop /3-functions at high temperature. Thus, in 
particular, the quartic coupling of (/?^ theory is shown to increase with 
temperature, in contrast to the behavior obtained in some previous 
studies. 
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1 Introduction 



Intuitively, one expects symmetries that are spontaneously broken at T = 
to be restored at high temperature. Examples of this behavior in condensed 
matter systems abound, and it was suggested many years ago [|[] that the 
same thing might happen in relativistic field theories. This idea was quickly 
backed up by calculations [0, ||, which essentially showed that for a general 
gauge theory with scalars and fermions the quadratic part of the effective 
potential goes like some combination of coupling constants times T^. For 
theories with simple scalar sectors this combination of couplings is quite gen- 
erally positive-definite, in which case the minimum of the effective potential 
lies at the origin of field space for sufficiently high temperatures. At the 
same time it was also pointed out however, that for certain models the 
coefficient of the term may be negative. In this case symmetries are not 
restored at high temperature, and in fact theories which are symmetric at 
T = may be broken for large T. While such theories are not totally generic, 
they are by no means contrived: all that is necessary is a sufficiently compli- 
cated scalar sector (two multiplets will often suffice) and some constraints on 
the relative sizes of the couplings in the theory. The later are usually weak 
enough so that the allowed values for the couplings occupy a sizable por- 
tion of parameter space. To be more precise: a necessary condition for high 
temperature symmetry breaking is the existence of negative scalar couplings 
in the Lagrangian. Such couplings, however, are automatically allowed in 
any model with at least two scalar multiplets, and consequently in almost 
all extensions of the standard model. The effect of these couplings on the 
phase structure of a theory depends on the details of the model and must be 
studied on a case by case basis. 

It was quickly realized that theories with symmetry non-restoration or 
multiple high temperature phase transitions have many attractive phenomeno- 
logical features. A partial list of applications includes the strong CP problem 
and GUT scale baryogenesis Q], the monopole problem ^, baryogenesis 
and dark matter [0, and inflation p. 

In this paper we will not discuss any particular physical application of 
models that exhibit an inverse phase structure. Instead, we will focus on the 
phenomenon itself, speciflcally on the validity of the one-loop calculation that 
predicts it in the flrst place. There are several reasons for doing this, perhaps 
the most important being that there seems to be a widespread suspicion that 
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symmetry non-restoration is merely an artifact of perturbation theory and 
not a true physical effect. While this is partially due to the counter intuitive 
nature of the phenomenon, there is some quantitative evidence to back up 
this claim. For example, a popular theory whose effective potential predicts 
an inverse phase structure in the one- loop approximation is the 0{N) x 0{N) 
model . However, subsequent calculations based on large N expansions and 
Gaussian effective potential techniques seemed to show that the symmetry is 
in fact restored at high temperature 0] . While the validity of these results is 
not clear to us, the fact that they incorporate some non-perturbative physics 
clearly raises the specter that the conclusions drawn from the loop expansion 
are erroneous. 

There are additional reasons to be concerned about the validity of the 
perturbative calculation of the effective potential. For example, consider the 
case were the theory is symmetric at T = 0, in which case perturbation 
theory predicts that the symmetry will be broken at high temperature. It 
is well known that the loop-expansion breaks down in the vicinity of the 
phase transition due to severe infrared divergences. Now we expect these 
infrared effects to strongly renormalize the coupling constants in the theory, 
which is important because the existence of the symmetry breaking transi- 
tion depends on a certain relationship between these couplings. It is clearly 
not inconceivable that the renormalization will destroy this relationship and 
circumvent the transition all together. Even in the case where one starts in 
the broken phase the situation is not entirely clear-cut. For while the loop- 
expansion predicts that this symmetry will remain broken at arbitrarily high 
temperature, it fails to include the effect of the temperature on the coupling 
constants. Again it is possible that the effective coupling constants appro- 
priate at high temperature do not obey the inequalities required to keep the 
symmetry broken. 

In light of all this it would be nice to study theories which naively admit 
an inverse phase structure using techniques that: one, capture enough non- 
perturbative physics to correctly handle the infrared problems in the vicinity 
of second order (or weakly first order) phase transitions; and two, unambigu- 
ously take into account the effects of temperature on the parameters in the 
theory. Fortunately such a technique has become available in recent years. It 
is based on an exact renormalization group (RG) equation for the Wilsonian 
effective action, which has the nice feature that it admits relatively simple 
approximations that nonetheless capture a good deal of the non-perturbative 
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physics[[ Using this tool we will study a Z2 x Z2 symmetric scalar field theory, 
the simplest model that exhibits an inverse phase structure in the one-loop 
approximation. As explained above, our goal will be to test the validity of 
the perturbative calculation and to establish if high temperature symmetry 
breaking is in fact possible. 

The paper is organized as follows. In section ^ we present a brief deriva- 
tion of the exact RG equation for the effective action and discuss the ap- 
proximate version that we will solve. To illustrate the techniques involved, 
section |] treats the high temperature phase transition of simple ^p^ theory. In 
section ^ we then focus on inverse symmetry breaking in the Z2 x Z2 model. 
Section ^ discusses the high temperature behavior of coupling constants and 
its relation to phase transitions. Our results are summarized in section ^ 



2 Flow Equation 

For completeness we present a brief derivation of the exact renormalization 



group equation. For details the reader is referred to the literature |]10|, |TT|, |12 



We work in D Euclidean dimensions and for simplicity consider a single real 
scalar field For a theory with action SYp\ a scale dependent partition 
function is defined as (dots represent contractions in function space) 

expIVA[J] =iV^P(^exp{-^(^- Ax'-<^-5aoM + J-<^}, (1) 

where Aq is an ultraviolet cutoff, is a J independent normalization factor, 
and 

is a free massless propagator times the ratio of two smooth (everywhere 
positive) cutoff functions Q^. This function has the properties 

6',(g, A) ^1 for g > A e 

and 

6^{q,K) ^0 for g < A - e. 



^See section S for details and references. 
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For later use we note that the sharp cutoff hmit is given by the Heaviside 
function: 

\ime,{q,A) = e{q- A). (3) 



The effect of the extra term in the exponential of (|l|) is to suppress the 
propagation of long wavelength modes (g^ <^ A^) while leaving the ultraviolet 
modes unaffected. In fact it is obvious from (|l]) that 



]iiaWA[J] = W[J] 



(4) 



the usual generating functional of connected Green's functions. The flow 
equation for Wa is obtained by differentiating ([^): 



OWaIJ] 
dA 



1 
2 



6Wa dA] 



6Wx 



5J dA 5J 
^ fdAlW'WAX 

IT 



V dA SJ5Jjj 

Now define a scale dependent Legendre effective action via 

1 



-ip-A],^-ip- Wa[J] + J-^ 



(5) 



(6) 



where 



5Wa 



In terms of Fa, (|^) becomes 
SFaM 



dA 

It is clear from (Hf) and 



1 

— tr 
2 



J_dAA 
Aa dA 



1 + A, 



S'Ta 



that 
1^0 



(7) 



(8) 



the generating functional of one particle irreducible (IPI) diagrams. In the 
opposite limit, A ^ oo, A^^ diverges, so that in (|l|) the "classical" approxi- 
mation to PFa[J] becomes exact. Thus 



Wa[J] ^ -S[^*] - -if* ■Al'-if* + J 



(9) 
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where 



5ip 

This says simply that as A ^ cxd, Wp^ becomes the Legendre transform of 
S + y-^l^ ■ if. But from (|) Ta + iy; ■ A^' ■ Lp is the Legendre transform 
of VFa, and since S" + ■ A^"*^ ■ is always convex in the limit A — oo, we 



obtain 



iimrAM = 5[v9]. (10) 

A — >oo 



We also have the approximate relation 

v^M = SkM. (11) 

which follows from (|TU]) provided the UV cutoff Aq is sufficiently large com- 
pared to all mass scales in S'aq • From (H) and (p!OD we see that Fa interpolates 
between the classical action at Aq and the effective action at A = 0. In fact, 
Fa has an interpretation as a Wilsonian quantum effective action obtained 
by integrating out purely quantum modes with momenta g > A ||T0| . 

Equation (|^) is exact, but much too complicated to be solved exactly. 
Its usefulness therefore hinges on the existence of sensible approximation 
schemes. We begin by writing 



f/A(^) + ^ZA(v^)(9V)' 

+Yt.{f,{d^ff){d^^f\. (12) 



Plugging this into (0), setting (p equal to a constant, and neglecting the ip 
dependence of we arrive at an approximate evolution equation for the 
effective potential: 

dV^{ip) 1 /• d^g a0.(g,A) 



dk 2 J (27r)^ OA 

X {e,{q, A) + e',{q, A)[Za - 1 + U'l{ip)/q']y' 

(13) 

(primes denote differentiation with respect to ip). Taking the sharp cutoff 
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limit and dropping an infinite field independent term we finally obtain^ 



dA 



K 



^A^-Mn 



A2 



(14) 



where Kc = Sd-i/ (2vr)^, and So is the surface area of a unit D-sphere. 

This is the equation we will study in the present paper, but before pro- 
ceeding a few remarks are in order. First note that ([T^) and ([llD have been 
obtained from the exact result (|^ by a truncation of the operator basis. 
While this is an uncontrolled approximation, the resulting equations have 
been successfully apphed to the study of phase transitions in two, three, and 
four dimensions, at zero and finite temperature, including the determination 
of accurate critical exponents 0, 0, |15|, |16|, [l^. They thus capture a large 
amount of the relevant physics. It is also possible to show that ( |T^ ) is simply 
the first term in a systematic expansion of the exact result [jl8|, which pro- 
vides some theoretical justification for its remarkable success. Finally, it is 
interesting to note that if one fixes Ua = Ua^ and Za = 1 on the RHS of (Q), 
and then integrates, the result is simply the well known one-loop effective 
potential. We emphasize, however, that there are no loop corrections to ([T^): 
the use of Ua on the RHS corresponds to a resummation of an infinite subset 
of diagrams, to all orders in the loop expansion. 

In this paper we will study the phase structure of four dimensional field 
theories at finite temperature. We will work with the sharp cutoff equation 
(p!4D, which, for the purpose at hand, has at least two advantages over its 
smooth cutoff counterpart: the flow equations are simpler, and there is no 
dependence on the particular choice of cutoff function. Using the sharp cutoff 
will also allow us to compare the two approaches since finite temperature 

studied using a si 
4 and setting Za 



phase transitions have been studied using a smooth cutoff in [|T5 
Specializing now to D 
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1 0, we are faced with the 



^Taking the sharp cutoff limit requires some care. We use the relation|10[ 



dA 



f{0e{q,A),A) 



5{A-q) / dtfit,q) 







(here f{9^,A) must be continuous at A = g in dependence on its second argument). 

reasonable approximation for the theories at hand, since the anomalous dimensions 
are known to be small . 
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task of including finite temperature effects. We begin by writing (jM]) as 



dK 



(2^ 



5{K-q) Hq^ + Ul{^)] 



The imaginary time formalism instructs us to replace 



(2vr) 



(2vr) 



and 



2 2,2 2,2 

q =g4 + q ^UJ^ + q_ 



where a;„ = 27mT. This yields 

dUA{f,T) 



dA 



2nT'Ag(^^]\n[A' + U';{v,T)], 



where 
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[27:t] 

n=0 



— )■ 

27iTj 



A 
2^ 



(15) 



(16) 



(17) 



and [x] is the greatest integer < x. We note that g{x) is continuous and 
piecewise differentiable (see Fig. p. 
It is easy to see that 



g(x) — as x — s> oo 
Stt 



so that 



reduces to 



9[/a(<^,T = 0) ir4,3 



OA 



A'' In 



A' + f/A(^) 



(18) 



(19) 



in the limit T ^ 0. Comparison with (|Tj) shows that this is indeed the cor- 
rect zero temperature equation. It is also instructive to look at the opposite 
limit, T ^ A. Using g{x) = x/27c for x < 2tt we obtain 



dA 



—TA' In 
2 



for T > A. 



(20) 



7 




5 10 15 20 

X 

Figure 1: The function g(x). 



On comparing with (11^ ) we see that this is simply T times the zero tem- 
perature flow equation for a three dimensional theory. Hence we already see 
evidence for dimensional reduction at high temperature. 

For later use we record here the analog of (|1^) for theories involving 
multiple scalar fields. The generalization is entirely straight forward, the 
result being 



9f/A(<^,T) K: 



dA 



^2'nT^Ag{^^ Indet 



A2 + 



3 J 



(21) 



where g{A/T) is given by (|T7|) as before. 

Before concluding this section we point out that there exists in the lit- 
erature an alternative method of including finite temperature effects in the 
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sharp cutoff case |]T9[. It amounts to replacing g by |q| in the delta function 



of ([151) . The result is 
OA 



oc Y: H^l + A' + U';{^,T)]. (22) 



The sum over Matsubara frequencies may now be done analytically, and 
the resulting equation involves no dependence on n, in contrast to ([T6|) . 
However, by cutting off only the three momenta and summing over all n 
one is neglecting the effect of integrating out the Matsubara modes with 
Un > A. on the effective potential. As noted following (H), it is precisely this 



kind of feedback that leads to the incorporation of higher loop effects. The 
equation derived from (^) has the additional unpleasant feature that it does 
not reduce to the zero temperature equation ([T^) in the limit T ^ pO[]. 



3 A(/?^ Theory 



In this section we will discuss briefly the finite temperature phase transi- 
tion for a simple Z2 symmetric scalar field theory. This model has already 
been analyzed using the Wilson RG approach [l^, and it is included 
here mainly to illustrate the method. It is worth noting, however, that the 
previous treatments differ from ours in the details of the implementation: 
13] uses a smooth cutoff, and in PO] a sharp cutoff is used only for the 



three-momentum, as discussed in the last paragraph of the previous section. 
Consider then a Z2 symmetric theory with Lagrangian 



C 



(23) 



where p = ^(p'^. Our RG equation for the effective potential (|TB|) becomes 



dUA{p,T) 
dA 



^ dUj,{p,T) ^ ^^d'UA{p,T) 



dp 



dp^ 



(24) 



We wish to solve (p^ ) subject to the boundary condition]^ U\g{p,T) = V{p). 
This is still rather difficult, and to make progress we expand Ua in a Taylor 

^See (|ll|). Note that consistency requires that we keep T <C Aq at aU times. 
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series about its minimum. This leads to an infinite set of coupled nonlin- 
ear ordinary differential equations, which may be approximately solved by 
truncation. 

Let us start with the symmetric regime, where the minimum of [/a is at 
the origin. We parameterize Uf^ in terms of its successive derivatives: 



m^(A,r) = f/P(p = o,r) 



^2„(A,T) = i-f/i")(p = 0,T) for n>2 



(25) 
(26) 



where U 



(n) 



A — 

derived by differentiating 
obtained in this way are: 



^"^"^ Evolution equations for these parameters are easily 



I) with respect to p. The first three equations 



dm? 
dXe, 



dA 



' 27r 

AT^ 

' 2ti 



6A4 



A2 + 



15Afi 



18 A/ 



A2 + w? 
28 As 



(A2 + n?f 
9OA6A4 



(27) 
(28) 



A2 + m2 (A2 + m?y 



+ 



72 A. 



(A2 + m?Y 



(29) 



Several remarks are in order at this point. First, note that the above 
expressions have an obvious interpretation in terms of one-loop Feynman 
diagrams. For example, (^81) says that the four point function receives con- 
tributions from the diagrams in Fig. ^. Figure ^(a) corresponds to the second 





(a) (b) 

Figure 2: Contributions to the running of A4. 
term in square brackets of (p3), which involves two 4-point couplings and two 
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"propagators" , while Fig. ^(b) corresponds to the first term in square brack- 
ets, involving the 6-point coupling and one "propagator". The other flow 
equations can be given similar interpretations. We also see how integrating 
out the high frequency modes generates new local interactions: even if one 
starts with a ip'' potential at Aq, higher order terms are created as the cutoff 
is lowered. Finally note that the equations describe naturally how particles 
with mass larger than A "decouple", in the sense that their effect on the 
running parameters becomes very small . 

Now let us derive the flow equations in the broken regime, where the 
minimum of the potential is at some po > 0. Here we parameterize [/a in 
terms of the location of its minimum and consecutive derivatives at that 
point. The flow equation for po{A,T) is obtained from the condition 



dp 



(po,T) = 



(30) 



Taking the A derivative of (|30[) we obtain (primes denote derivatives with 
respect to p): 



dU: 



dk ^ dA 



so that 



dpo 
dA 



dU^J_ 
OA U'i 



(31) 



(32) 



(all ?7a derivatives are evaluated at po, of course). The other parameters are 
defined as in (PB|): 

(33) 



A2„(A,r) = ^f/f (po,T) 



for n>2. The flow equations for these parameters are obtained by differen- 
tiating ( p^ ) with respect to p and evaluating the result at po- The first three 
are: 



(iA4 



AT' 

' 2n 



dpo 
dA 

A 



AT' 



A 



6 A4 + 12 poAe 
A2 + 4poA4 



1 
2AI 



'[t 



27T 

15A6 + 24poA8 

A2 + 4poA4 2(A2 + 4poA4)' 



(6A4 + 12poAe 



\2n 



(34) 



.3A.^ . (35) 
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We begin our study of the high temperature phase transition by com- 
puting the effective potential in the simplest possible approximation. This 
amounts to parameterizing Ua in terms of its first two derivatives only, 
i.e., we set A„ = for n > 6. Our strategy is then as follows: we de- 
fine the theory by specifying "renormalized" zero temperature parameters 
at A = T = 0. We may choose the zero temperature theory to be in the 
symmetric regime by picking mf. = m^(0,0) > or in the broken regime 
by picking por = Po(0,0) > 0. Depending on the chosen phase, we then 
use either , (pS]) or (P^,(P5|) to numerically integrate up in A to some 
Aq ^ rrir or y^por- This "integrating up" is done at T = 0, and it serves 
merely to provide us with the "bare" parameters that define the action Sa^. 
We note that if or po go to zero at some intermediate A, then we con- 
tinue the integration with the set of equations appropriate to the new phase, 
taking the values of the parameters at A as initial conditions. 

With the parameters at Aq in hand we are now ready to study the ef- 
fects of finite temperature. This is done by fixing T 7^ in the evolution 
equations and running down from Aq to A = 0. Repeating this for differ- 
ent values of T allows one to determine the temperature dependence of the 
renormalized parameters. Figures ^ and ^ show the evolution in A of the 
parameters at various temperatures. The initial conditions are such that the 
zero temperature theory is in the broken phase. From Fig. |^ we see that po is 
quadratically renormalized for large A and approaches a constant as A — >^ 0. 
As the temperature is increased, the asymptotic value of po decreases, until 
at Tc we have po(A, Tc) as A 0. Above po goes to zero for finite 
A, and we continue the evolution in the symmetric regime. Figure ^ shows 
the corresponding behavior of A4. At T = we observe the expected log- 
arithmic evolution. Above and below Tc, where all modes are massive, the 
running of A4 is effectively stopped when A becomes much smaller than the 
relevant mass scale. At Tc, however, A4 runs to zero. This is, in fact, ex- 
pected, because at the critical temperature the theory is scale invariant and 
the phase transition is second order. As a consequence the parameters run 
according to their canonical dimension, which is one for A4 at small A since 



28 As + 40 pqAio (15 Ag + 24 ppAg) (6 A4 + 12 ppAe 
A2 + 4poA4 



(A2 + 4poA4) 



+ 



(6A4 + I2P0A6: 
3 (A2 + 4P0A4) 



(36) 



12 



3 



2.5 



1.5 



0.5 



Po(T=0) 
Po(T=2) 
Po(T=2.92) 
Po(T=4) 
m^(T=4) 



0.001 



0.01 



0.1 
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Figure 3: The evolution of po at various temperatures. For T > T^, the 
running of the mass term in the symmetric regime is also shown. 



to the long wavelength modes the theory looks effectively three dimensional 
at high temperature^ It is the fact that the Wilson RG approach correctly 
takes into account the strong renormalization of the coupling in the vicinity 
of the phase transition that allows one to accurately describe the physics in 
this region, including critical exponents |]T3|, |T^. This is to be contrasted 
with the loop expansion, where higher order terms go like \i{r / M. to some 
power, with M. the appropriate infrared cutoff. Since ^ at Tc, the 
expansion necessarily breaks down near the phase transition. 

Figure ^ shows the temperature dependence of the "renormalized" (i.e. 
A = 0) parameters. The critical temperature is given by T^jp^r ~ 8.60, 
whereas naive one- loop perturbation theory predicts T"^ / por = 8 p|, H. In 
light of the fact that the loop expansion breaks down near the phase tran- 
sition this good agreement may seem surprising, but it is in fact expected: 



^This is explained in more detail in the appendix. 
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Figure 4: The evolution of A4 at various temperatures. 



since perturbation theory is vahd both far above and far belowQ the critical 
temperature, one can show that the error in the naive value of Tc is of order 
XiTc 0. The above results may also be compared with ||T^, where finite 
temperature ip'^ theory was treated using a smooth cutoff flow equation. 

Before going on to theories involving multiple scalar fields a discussion of 
the validity of our approximation method is in order. We will not comment 
on the initial step of obtaining (|14[) from (|^) by discarding the momentum 



dependence; for this we refer the reader to the literature [|T0|, |18[. What 
we do want to address is the effect of truncating the infinite set of coupled 
differential equations. Instead of discussing the formal aspects of the problem 
(see for example [Tl|, |T^), we will focus on a simple practical test: if the 



truncation is to be sensible at all, then the effect of including additional 
terms must be, in some sense, "small". 

To see if this is the case we derived flow equations analogous to 



^As long as there are no Goldstone modes. 
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Figure 5: The temperature dependence of the renormahzed parameters. 



and (|34 )~(p6D up to and including A14. Solving this set of seven coupled 
equations requires a change in tactics. Recall that before we were able to 
specify the renormahzed (A = 0) parameters and then integrate up to ob- 
tain the "bare" parameters at Aq. This is no longer possible, because the 
newly included couplings A6-A14 correspond to irrelevant operators. This 
means that their values in the infrared are fixed in terms of the values of the 
relevant parameters and A4, and consequently they may not be chosen 
independently |]22[| . 



To isolate the effect of including the additional terms we therefore proceed 
as follows. Start at Aq and set A6-A14 arbitrarily to zero|. Then fine tune 
Po(Ao) and A4(Ao) until flowing to the infrared produces por = Po(A = 0) = 1 
and X^r = A4(A = 0) = 0.1. At this point we know we are dealing with 
the same theory as before, and we may now switch on the temperature and 
find Tc. This procedure was carried out for several different values of A4, 

^Any other "natural" value would do just as well, where "natural" means 0(1) divided 
by enough powers of Aq to get the dimension right. 
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Table 1: T^^/por ioi various approximation methods and coupling values. 





Calculation Method 


0.01 


0.1 


1.0 


1-loop Pert. Th. 


8 


8 


8 


Wilson RG up to A4 


8.29 


8.60 


10.37 


Wilson RG up to A 14 


8.28 


8.55 


9.69 



the results being summarized in table |l|. We see that at least for moderate 
couplings the effect of including parameters beyond A4 is small. This is in 
accord with the expectation that as long as the effective potential is not too 
steep the first two terms in the Taylor expansion should accurately describe 
its shape around the minimum. 

4 Inverse Symmetry Breaking 

The simplest model that exhibits an inverse phase structure is a Z2 x Z2 
symmetric scalar field theory described by the Lagrangian 

c = ^{d,vr+l{d,xr-v{v,x), (37) 

where 

Vif, X) = ^m If +2f^X + -f'P + -fx - X ■ (38) 

Note that, if the A's are all positive, we require 

W > Kx (39) 

for boundedness. If one calculates the finite temperature effective potential 
for this theory to one loop and then expands the result about the origin one 
finds that the quadratic part at high temperature is 

+ ^(-"' + ^(3A,-A,,)]x^ (40) 
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It is easy to see that this form allows a myriad of different symmetry breaking 
patterns. Depending on the relative size of the couplings and the signs of 

and /i^, one can have symmetry breaking, restoration, or non-restoration 
at high temperature. Although less obvious, there are even more exotic 
possibilities: for example, one can have a symmetry broken or restored only 
for an intermediate range of temperatures ||23|| . 

Let us now analyze this model using the Wilson RG formalism developed 
in the previous section. Because of the Z2 x Z2 invariance the Lagrangian 
may be written as 



where 



and 



c 



(41) 



(42) 



(43) 



We proceed as in the previous section. The RG equation ( pT] ) for this model 
is 



r-l±27iT'Kg(^^] IndetA^, 



where 



M 



92 {/ 



dpdC 



dpdC 

^ dc, ^ ac^ 



(44) 



(45) 



The potential [/a is again parameterized by its first two derivatives at the 
minimum, but this time we must distinguish four cases: the minimum may 
lie at the origin, on either axis, or between the axes. In all four cases the 
couplings are defined via 



A^(A,T) 
Ax(A,T) 



2 9p2 
2 dC^ 
2 dpdC 



(46) 
(47) 
(48) 
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The other two parameters used to describe the potential depend on the lo- 
cation of the minimum. If it is at the origin, we simply define 



dUA{0,0,T) 
dp 

dUA{0,0,T) 



(49) 
(50) 



If the minimum is at po on the p axis, we use 

dUA{po,0,T) 



/(A,r) 



and Po to parameterize Ua. The flow equation for po is derived as in 
Similarly, if the minimum is on the ( axis. 



(51) 



dUAiOXo,T) 
dp 



(52) 



and Co are used. The final possibility is that the minimum lies between the 
axes. In this case we parameterize in terms of its location in field space, 
(po,Co)- The flow equations for these coordinates are obtained by taking 
total A derivatives of the two expressions 



and 



^^a(Po,Co) 
dp 

dUA{po,Co) 
dC 







. 



(53) 



(54) 



This yields two simultaneous algebraic equations for dpo/dA and d(o/dA, 
which may be solved to give 



dpo 
dA 

dA 



\ ± ( 9Ua\ 
1 -^X dp \ dA J 



(po,Co) 



, X d_ { duA 

'^'fix ac \ dA J 



(po.Co) 



I'^f ac \ dA J 



-^ X ± 

(po,Co) ^ 9P 



(dUA\ 

{ aA J 



(po,Co) 



X X — X2 



(55) 



(56) 
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As before, the right hand sides of the evolution equations are expressed 



in terms of our parameters by differentiating (44) with respect to the fields, 
evaluating the result at the appropriate minimum, and then dropping all 
terms with three or more derivatives. For example, the flow equations for 
the symmetric regime (i.e. minimum at the origin) are as follows: 



dm? 

"dK 
dil_ 

dA 

dXtp 

UK 

dA 
dA 



AT^ /A 
A 



2n 
AT2 



6A„ 



2A 



2n 



A^ + m? 
6Av 



A2 + |i2 

2 A^^ 



A2 + ^2 A2 + m? 



Ar2 




^A 


27r 


'A 


vT 


AT2 




^A 


27r 




vT 


Ar2 




^A 


27r 




vT 



18 A^ 



(A2+m2)2 (A2 + ;i2)2 



18 A^ 



+ 



(A2+/i2)2 (A2 + m2)2 



(57) 
(58) 
(59) 
(60) 



6 A<^A<^^ 



+ 



6 A^A^^ 



(A2+m2)2 (A2 + ;i2)2 



'Kx 



(A2+m2)(A2 + /i2) 



(61) 



There are similar, albeit more complicated, sets of equations for the other 
three positions of the minimum. 

We begin our numerical study of the model by investigating the phe- 
nomenon of inverse symmetry breaking. By this we mean choosing the pa- 
rameters so that the vacuum is symmetric at T = but asymmetric at high 
temperature. From (EOl) it is clear that according to the one-loop result this 



happens if, say, A^,^ > 3A(^. The critical temperature is predicted to be 
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^ipx 3A(^ 



(62) 



independent of A^^ and /i^. 

The flow equations are solved as in the previous section. Renormalized 
parameters are chosen at T = A = in such a way that the vacuum is sym- 
metric. The equations are then integrated to some Aq much larger than all 



(p9) then requires that > X^p^^ so the other symmetry cannot be broken at high 



T. 
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physical masses. Since the minimum moves away from the origin during this 
process it was necessary to develop an algorithm that automatically switches 
to the correct set of equations depending on the location of the minimum. 
Once the parameters at Aq are known the temperature is switched on and 
the equations are integrated back down to A = 0, yielding the renormalized 
finite temperature values for the masses, vacuum expectation values, and 
couplings. 




Figure 6: The masses and por as a function of temperature. 

Figures ^ and ^ show the behavior of the renormalized parameters for a 
typical case. The parameters are chosen so that naively (i.e. according to 
(^OD) the Z2 symmetry of the (/9-field is broken at high temperature. From 
Fig. P we see that this picture is confirmed by our RG approach. As the 
temperature increases m^(T) decreases and eventually hits zero at Tc ~ 49. 
Both ii^(T) and por(T) are proportional to at high temperature. Fig- 
ure ^ shows the corresponding behavior of the couplings. As in the if^ case 
discussed in the previous section, there is dramatic renormalization in the 
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vicinity of the phase transition^. 
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Figure 7: The couphngs function of temperature. 

Figure ^ compares the numerical value of the critical temperature with the 
naive prediction, equation (^), for several values of the cross-coupling. Even 
for the relatively small couplings chosen there are significant differences for 
all values of A^^. The best agreement occurs for large values of A^^, but even 
here the difference is about 15%. (Note that (^) requires that A^^ < 0.0173. 
Numerically the system becomes unstable shghtly later, at A^^ ~ 0.018.) As 
A^^ decreases, the deviation between the naive prediction and the numerical 
result increases rapidly. For A^^ = 0.0105, the difference is more than a 
factor two. For AJ^^ = 0.01, perturbation theory predicts Tc/rrir ~ 109. The 
numerical solution shows that there is no phase transition at all for these 
values of the parameters, and that m^(T) is in fact an increasing function of 
T at high temperature. 

That the RG calculation gives a significantly higher critical temperature 

^"^The behavior of the couphngs near is discussed in more detail in the appendix. 
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Figure 8: Comparing the critical temperature obtained from one-loop per- 
turbation theory with that obtained by integrating the RG equations. (AJ^ = 
0.003, = 0.1, ^4 = 1) 

than perturbation theory is due to a combination of two factors. First of all, 
the renormalized couplings decrease with temperature, which raises T^. This 
effect, however, is also present in the simple (f'^ case, and as was shown in 
section 0, there the two methods agree quite well. It alone therefore cannot 
account for the difference in critical temperatures. In fact, the dominant 
reason for the discrepancy is most easily seen in Eq. (p?]). The "propagator" 
in the second term is -|- /x^, and /i^ is a rapidly increasing function of 
temperature. Hence at high T the second term is suppressed compared to 
the first, which significantly diminishes the effect of the cross-coupling A^^, 
and consequently raises Tc. 

We have just seen that the temperature dependence of /i^ plays an impor- 
tant role in determining Tc. Since the leading correction to /i^ is proportional 
to we may therefore suspect that the critical temperature also depends on 
the value of this coupling, even though the one-loop result ( |62|) predicts that 
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it does not. The question is settled by Fig. ^, which shows the dependence 
of Tc on A;^. The effect is clearly quite large and increases rapidly as the 
coupling increases. Note that > 0.0521 is required for boundedness. 

150 
140 
130 
120 
110 
100 

90 

80 

70 

60 

50 

Figure 9: The dependence of the critical temperature on A^. (A^ = 0.003, 
X:^^ = 0.0125, /i2 = 1) 

Figures H and ^ indicate that the critical temperatures obtained from the 
RG approach are consistently higher than those gotten from the loop expan- 
sion, and that certain sets of couplings don't produce a phase transition at all 
even though X^^^ > SXip. This may lead one to the conclusion that the region 
of parameter space which yields an inverse phase structure is significantly 
smaller than that predicted by the one-loop result. We emphasize that this 
is not the case. The reason is again the presence of the "propagators" in the 
flow equations. As explained above, for the case at hand these factors reduce 
the significance of the X^p^ term, which raises Tc and suppresses the phase 
transition. But at the same time the suppression of A^^^ tends to stabilize 
the system, which means that there is additional parameter space beyond 
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Table 2: T(./mr for various approximation methods and couplings (/x^ ^ 1 in 
all cases). "Up to dim. n" means that operators of dimension higher than n 
were discarded in the flow equations. 



K = 
f- = 


2.78 X 10-3 
1.00 X 10-^ 
1.25 X 10-2 


1 X 10-2 
3 X 10-1 
5 X 10-2 


1 X 10-4 
3 X 10-3 
5 X 10-4 


1-loop Pert. Theory 


53.7 


24.5 


245 


RG up to dim. 4 ops. 


67.9 


33.6 


257 


RG up to dim. 8 ops. 


67.9 


33.5 


257 



the limit ( P^ ) for which transitions occur. The net effect is then a shift of 
the inverse-breaking region rather than a major reduction. 

To make sure that the numerical results do not depend significantly on 
the order of truncation of the evolution equations it is necessary to study 
the effect of including higher dimensional operators. To this end we enlarged 
our system of equations to include the couplings corresponding to operators 
of dimension less than or equal to eight, e.g., ^9^, (^"^x^i X^i stc. This results 
in a system of sixteen coupled first order nonlinear differential equations. As 
discussed at the end of section |^, the addition of irrelevant operators causes 
the complication that one can no longer simply choose the renormalized 
parameters. Instead, one must fine-tune the bare parameters at Aq to achieve 
the desired values for the masses and relevant couplings in the infrared. Aside 
from this nuisance the solution of the equations proceeds as before. The 
results are summarized in table |^, which shows that the effect of the added 
terms is very small. We conclude that discarding operators of mass dimension 
higher than four is a valid approximation for the purpose at hand. Table ^ 
also shows that for very small values of the couplings the numerical result is 
in reasonable agreement with the one-loop prediction, at least if \^p^ is safely 
larger than 3A(^. For the couplings in column 2 the two methods differ by 
almost 40%, while the discrepancy is down to 5% for the values in column 
three. 
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5 The Coupling Constants at High Temper- 
ature 



In this section we discuss the temperature dependence of the couphng con- 
stants in the high T hmit. We begin with the simplest case, that of a single 
scalar field. 



5.1 Xip^ Theory 

The high temperature behavior of the quartic coupling A has been the subject 
of controversy for some time. Several authors have found A decreasing with 
T others find A approaching a constant as T ^ oo and still others 



find that A is increasing at high temperature . 

In order to shed some light on the issue we integrated the RG equations 
discussed in section ^ for a wide range of temperatures. The results are shown 
in Fig. |10|. We see that the coupling constant rapidly decreases at first, but 
eventually turns around. From T f» 40 onwards the dependence is clearly 
approximately logarithmic. We note that at least qualitatively this behavior 
agrees with that obtained in . 

In the literature one frequently encounters the notion that the effect of 
high temperature on coupling constants may be incorporated by "running" 
the couplings according to their one-loop /3-functions using the temperature 



as the scale p8[. In our case this amounts to assuming 

Al , (63) 



dX/^j- 9 2 



dt 87r2 
where t = ln(T/To), so that 

^"^^^ = 1 - 3f.Al"Sln(T/To) • 

The justification for this procedure is usually based on rather vague argu- 
ments concerning the average momentum transfer during collisions of parti- 
cles in the heat bath being 0{T), so that a coupling constant at that scale 
ought to be appropriate. While this is certainly not unreasonable it is hardly 
a convincing argument, and as noted above a more detailed analysis has 
produced many other types of behavior. To see how well Eq. (Bl) describes 
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Figure 10: Temperature dependence of the quartic coupling constant. Note 
the logarithmic scale on the temperature axis. {X^r = 0.1) 

the high temperature evolution we have plotted Eq. together with the 
numerical integration in Fig. |TT]. The one-loop result is normalized so that 
the two methods agree at T/m^ = 100. It is clear from the figure that the 
high temperature dependence of A4r.(T) is described quite well by the one- 
loop /3-function. After running T over three orders of magnitude the results 
differ by about 1.5%. It is also clear, however, that the difference becomes 
significant if one runs over many orders of magnitude, for example to the 
GUT scale. In addition we see that using the one-loop equation is only rea- 
sonable if one knows the correct "initial condition" at high temperature, in 
our case XiriT) ^ 0.085 at T/rrir = 100. Starting the evolution using Eq. 
( |6^ with A4,.(T = 1) = 0.1 is obviously not satisfactory, and obtaining the 
proper initial condition requires an analysis such as the one presented in this 
paper. 

It is worthwhile noting that the general features of the high temperature 
behavior of the coupling can be understood by simply looking at its RG 
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Figure 11: The high temperature evolution of A4r(T) obtained from the RG 
compared to the behavior obtained by using T as the scale in the usual one- 
loop /3-function. The initial condition for the one-loop integration is chosen 
so that the two methods agree at T=100. (A4r = 0.1) 



equation, obtained in section ^. For Ag = and A ^ T and \m\, Eq. (^) 
reduces to Eq. (|63D. Thus the evolution for very large A is simply given by the 
one-loop /3-function. Now for T 3> rrir the propagator in (pSD contains a large 
thermal mass, which means that the running effectively stops for A < y/X^T. 
The evolution for A <^ T is thus strongly suppressed, and all that remains is 
the one-loop-like running for large A plus some threshold effect for A ~ T. 
At a higher temperature the thermal mass is larger, causing the running of 
A4 to freeze out at larger A. The contribution from the threshold region 
is approximately the same as before, and so the net effect of increasing T 
from Ti to T2 is essentially equivalent to integrating the one-loop /3-function 



between the two temperatures. This explains the correspondence in Fig. 11 



The slope of the numerical curve is slightly steeper because the quantum 
corrections to the mass go like — A4A^, which reduces the denominator in 
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and makes A4 flow somewhat faster for large A than predicted by Eq. 
30. 

Having just explained why X4r{T) should increase with temperature we 



find ourselves having to account for the initial decrease seen in Fig. |T0|. This 
decrease occurs for values of the temperature low enough so that the suppres- 
sion of the evolution due to the thermal mass is not yet important. Instead 
the running is dominated by the range A < T, where the prefactor on the 
RHS of Eq. (p^), AT^(yf(A/T), is a rapidly increasing function of T. Thus, for 
moderate values of T, the coupling runs faster at higher temperature, which 
explains the initial decrease in X^riT). 

We end this section with a few remarks. First of all, we have checked that 
the high temperature behavior of the coupling is not altered by the inclusion 
of higher dimensional operators in the flow equations. In fact, the effect of 
these operators decreases rapidly as one gets further away from the critical 



theory. Second, we point out that the author of using an "improved" 
loop expansion approach, has found the high temperature evolution of the 
coupling to be similar to ours, i.e., approximately as predicted by the one- 
loop /9-function with some temperature dependent correction. Finally, note 
that the calculations in this section were done for a theory that is symmetric 
at T = 0. This choice has no effect on the high temperature behavior of the 
coupling, for at very high T the theory does not care whether or not there 
was a phase transition at much lower temperature^^. For the temperature 
dependence of the coupling constant in an initially broken theory at moderate 
values of T see section ^. 

5.2 Z2 X Z2 Model 

We will now study the high temperature behavior of the coupling constants 
of the Z2 X Z2 model discussed in section |. There are two possible phases at 
high temperature: either the theory is completely symmetric, or one of the 
two symmetries is broken. We have checked that in both cases for generic 
values of the zero temperature parameters the evolution of the couplings at 
high temperature is approximately given by their one-loop /3-functions, viz., 

§ - gi^K + AjJ. (65) 



^This was also checked numerically. 
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dt 



d\ 



dt 



— [9Aj + AjJ 
' ■[3(A^ + AJA 



Stt^ 



4A 



(66) 
(67) 



where t = ln(T/To). The reason for this behavior is essentially the same as 
in the ip'^ case discussed above, except that the situation is complicated by 
the presence of multiple mass scales. An example of the evolution of the 




1000 



10000 



100000 



Figure 12: The high temperature evolution of the couplings. "RG" refers 
to the Wilson RG approach of this paper while "1-loop" refers to the nu- 
merical solution of equations (|65|)-(|67|). The initial conditions for the one- 
loop integration are chosen so that the two approaches agree at T = 500. 



{ml = = 1, a; = 0.1, A^^ = 0.003, A^^ = 0.016) 



couplings is shown in Fig. |T^, which also presents the running according to 
the one-loop /^-functions for comparison. We see that at least over the limited 
temperature range presented the two methods agree reasonably well. Note 



that the renormalized couplings in Fig. |12| are chosen so that the theory is 
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symmetric at T = but broken at high temperature. The transition to the 
broken phase occurs at Tc ~ 50. 

We conclude this section by discussing an interesting phenomenon related 
to the evolution of the couplings. Consider choosing renormalized couphngs 
which satisfy 

> > . (68) 

From (|65D-(|6^) it is clear that according to the one-loop /^-functions the cou- 
plings will flow in such a way that eventually X,pxi^) > 3A;^(t). For suitable 
initial values this occurs when all couplings are still small, so that the loop 
expansion ought to be valid. This means that the symmetry of the theory 
should be given by Eq. (^), and we see that the running of the couphngs 
can alter the symmetry of the theory at high temperature. This effect can be 
exploited for model building, in that it provides a "natural" way of having 
phase transitions at very high temperature. For example, one can have sym- 
metry restoration set in above the GUT scale without having to fine tune 
the couplings to one part in 10^^. Similarly, one can construct models that 
are broken at low temperature, get restored at some ~ 0{m/y/\), and 
then get rebroken at a much higher temperature. For a practical application 
of this kind of effect to the strong CP problem and the baryon asymmetry 
see i. 

The above discussion was based on the assumption that the couplings 
evolve at high temperature according to their one-loop /^-functions. We have 
already stated that this is approximately true, and we therefore expect that 
the evolution may indeed reverse the inequahty 3A-^ > X^^^ at high temper- 
ature. We wish to investigate the question if this really leads to symmetry 
breaking, as predicted by (|40|) . Unfortunately the errors inherent in the 
numerical solution of our evolution equations do not allow us to vary the 
temperature over the 15 or so orders of magnitude necessary to achieve a sig- 
nificant change in the couplings. We therefore focus on the situation where 
one starts with A^^ less than, but very close to, 3A^, so that even a moderate 
change in T can alter their relationship. 

Bearing in mind the above restrictions we choose A^ = 0.1, A^ = 1 x 10~^, 
and A^^ = 2.99 x 10~^. If one numerically integrates equations (|65D -(|6^) 
using these couplings as initial conditions at t = 0, the one- loop calculation 
( ^Of ) predicts a symmetry breaking phase transition at T ^ 450. The results 
of the Wilson RG integration are presented in Fig. |TB|. It is clear from 
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Figure 13: The temperature dependence of two combinations of parameters in 
the Wilson RG approach. uliT) is rapidly increasing long after the couplings 
satisfy the naive inequality which predicts symmetry breaking, (m^ = /i^ = 
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the figure that the coefficient of the term in /i^(T) is not proportional to 
3A^ — A^^ as predicted by Eq. (|40|) . In addition we see that even long after the 
perturbative coefficient has turned negative, uliT) is still a rapidly increasing 
function of temperature. The phase transition predicted by the one-loop 
result is absent. This is not to say, of course, that a phase transition may 
not take place at much higher temperature. After all /i^(T)/T^ is decreasing, 
and if this continues the symmetry will indeed by broken at very large T. 
As stated above, numerical errors prevent us from investigating this region 
directly. 



Finally, we point out that for the initial conditions of Fig. |T3|, Ar,(T) does 



not run according to its one-loop /^-function over the temperature range con- 
sidered. In fact, A^(T) is decreasing very slowly even at T = 10^. This shows 
that while generically the high temperature behavior of the couplings is given 
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by their one-loop /^-functions (Fig. this need not be so. The "anomalous' 
behavior of A^(T) for the case at hand is due to the presence of two vastly 
different mass scales at high temperature. When m^(T) and jJ^liT) differ by 
several orders of magnitude, some terms in the flow equations "decouple" 
much earlier than others. Under these circumstances our approach is not ex- 
pected to reproduce the one-loop result as the latter is blind to the presence 
of mass scales. 

6 Conclusions 

In this work we investigated the high temperature phase structure of a simple 
two scalar theory by solving approximately a non-perturbative RG equation 
for the effective potential. The solution consisted of assuming a polynomial 
expression for the potential and then numerically integrating the resulting 
coupled flow equations for the coefficients. Our main result was obtained 
in section ^: according to our non-perturbative method high temperature 
symmetry breaking (or non-restoration) does exist. In addition, we found 
that the phenomenon takes place roughly for those values of the couplings 
that satisfy the inequalities obtained from perturbation theory. The total 
volume of parameter space that yields high temperature symmetry breaking 
is only slightly reduced compared to the one-loop prediction. 

We also saw that the critical temperature obtained from perturbation 
theory does not agree particularly well with our numerical results, even for 
reasonably small couplings (< 0.1). The reason for this is that the negative 
cross-coupling between the fields which drives the symmetry breaking must 
by somewhat larger than the value predicted by the loop expansion. This 
leads to a situation where the perturbative estimate for the critical tem- 
perature becomes totally unreliable as the cross coupling gets close to the 
boundary region. If one stays far away from the boundary the perturbative 
and numerical results approach each other to within 15%, which is about 
what one would expect. We also demonstrated that the critical temperature 
depends significantly on the value of the quartic coupling of the unbroken 
field, in contrast to what is predicted by the one-loop calculation. 

In section || we discussed the behavior of the coupling constants at high 
temperature. In the tp^ case the quartic coupling was shown to increase 
at very high temperature, but not before being significantly reduced during 
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an intermediate regime. For very large T the evolution was shown to be 
approximately given by the one-loop /3-function of ordinary perturbation 
theory. 

For the Z2 x Z2 theory we focused on the effect of the running couplings 
on the symmetry of the theory. It was shown that generically the couplings 
evolve according to their one-loop /5-functions at high temperature. However, 
we also demonstrated that it is not correct to draw conclusions about the 
symmetry of the theory based on the perturbative formula for the thermal 
mass and the running of the couplings. Naively this idea can be used to 
produce "natural" very high temperature phase transitions induced by the 
evolution of the couplings. Our RG approach shows that this fails for at least 
three reasons. First of all, even for generic initial conditions, the couplings 
evolve according to their one-loop /3-functions only at very high T, and they 
are significantly renormalized by the time they reach this region. Because 
the couplings evolve very slowly after the initial decrease it is important to 
have these "initial values" in order to start the high temperature running. 
Using the zero temperature parameters would not suffice to get an order 
of magnitude estimate of Tc even if the rest of the reasoning was correct. 
Second, any attempt to cause high temperature symmetry breaking induced 
by the running of the couplings requires that one of the quartic couplings be 
significantly larger than the others. This introduces two widely different mass 
scales (the two thermal masses), with the result that some of the couplings 
don't evolve according to their one-loop /^-functions even at very large T. 
Lastly, and most importantly, we showed that even when the couplings have 
evolved to fulfill the naive inequality that predicts symmetry breaking, it 
does not happen. 

The above remarks are best illustrated by the example discussed in section 
1^. For the renormalized parameters chosen there the one-loop evolution of 
the couplings (Eqs. (|65|)-(|67D) combined with the perturbative formula for 
the thermal mass (Eq. (PP])) predicts a symmetry breaking phase transition 
at T ^ 450. In contrast. Fig. ^ shows that the relevant thermal mass is a 
rapidly increasing function of temperature even at T = 10^. It is possible 
that the symmetry will be broken at much higher T, but if this happens it 
will be at a temperature many orders of magnitude larger than that predicted 
by perturbation theory. 

Finally we would like to comment on the generality of our conclusions. All 
of our numerical results were obtained for the simple Z2 x Z2 model. However, 
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the differences between our RG approach and the standard one-loop treat- 
ment arise because our method correctly takes into account the decoupling 
of massive particles from the theory as the infrared cutoff is lowered. This 
feature is intrinsic to the method and independent of the particular model 
studied, which leads us to expect similar deviations from the perturbative 
results for other models. 
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Appendix: Critical Behavior 

The main emphasis of this work does not rely on the details of the phase 
transitions we have studied. Rather, we have been interested in establishing 
the existence of symmetry breaking phase transitions, and in showing that 
the broken symmetry state can persist at arbitrarily high temperature. In 
this context it is only of tangential interest as to what critical exponents 
characterize the transition or even whether it is first or second order. In 
fact, we could have studied the persistence of the broken state by starting 
in an asymmetric vacuum at T = 0, in which case there would have been no 
phase transition at alQ Similarly, the evolution of parameters at very high 
temperature does not depend on the detailed dynamics of a phase transition 
that has taken place at much lower T. The nature of the phase transition 
is important, however, in understanding the behavior of the renormalized 
parameters near Tc, shown in Figs. |^-|^. For this reason we include here a 
brief discussion of the phase transitions studied in sections ^ and ^. Detailed 
applications of Wilson type RG equations to critical phenomena can be found 

in ig,|n], [15 0,0. 

We begin with the well known case of Ay?^ theory. In order to study 
the character of a transition it is helpful to cast the flow equations into 
scale invariant form by using the proper dimensionless couplings. To achieve 
this consider the theory at high temperature. Once A < 27rT all non-zero 

^^So long as the couplings obey the correct inequahties, see the discussion above (|6^). 



34 



Matsubara modes have been integrated out and we are left with an effective 
three dimensional theory for the zero mode (see (pOD ). The coupling in this 
theory, which has dimension one, is A4T. The appropriate dimensionless 
mass and coupling parameters of the effective theory are hence 

.,A,r).!!%I) (69) 

and 

MA,T) = M^. (70) 

Rewriting equations (|27|) and (pS]) in terms of n and h and using t = ln(Ao/A) 
we obtain the desired form: 

dK ^ 3h 

IE = 

dh _ ^ 9/1^ 

H ~ ^"27r2(l + «:)2 ^^^^ 

(here we have truncated by setting An=0 for n > 6, as before). Equations 
(fflD and ( [72D have two fixed points: the Gaussian (trivial) fixed point at 
h-^ = hi, = 0, and the Wilson fixed point (WFP) at = —1/7 and = 
87r^/49. Linearizing around the fixed points one finds that the former is 
completely unstable in the IR {t 00), while the latter is a saddle point. 
Hence in order for a flow to end up at the WFP one needs to fine-tune one 
linear combination of UV couplings, which in our case amounts to choosing 
T = Tc. Consequently k ^ and — >■ /i^ as t — >■ 00 at the critical 
temperature, which means that the transition is second order and that 

m^(Tc) = hm m^(A, ~ lim = , (73) 

Ah 

X,r {Tc) = hm A4(A, TJ ~ ][im ^ = 0. (74) 



This explains the behavior observed in Figs. near Tc. 

We point out that the above conclusions are independent of the parame- 
terization used for the the flow equations, as they should be. For example, 
consider starting with the flow equations in the broken phase, (0) and (|53|). 
The appropriate dimensionless parameters of the effective three dimensional 
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theory are now k = po{A, T)/AT (recall that the field has dimension one-half 
in three dimensions) and h = A4(A,T)T/A. In terms of these variables (0) 
and ( P^D take a scale invariant form similar to ( [7TD and (|72p. One again 
finds two fixed points, the Gaussian at = 3/47r^, /i^ = and the WFP 
at = l/47r^, h^, = 27r^. The former is again completely unstable and the 
latter is again a saddle point, so our above conclusions remain valid. 

We now turn to the two scalar theory of section ^. Just as above we 
consider the effective three dimensional high temperature theory and rewrite 
the flow equations ([57|)-(|6TD in terms of the dimensionless variables = 
m2(A,T)/A2, = /i2(A,T)/A2, = A^(A,T)T/A, = A^(A,T)T/A, 
and htp^ = A^^(A,T)T/A. This results in the following set of scale invariant 
equations: 



~dt 
~dt 
dt 
~dt 



dh. 



dt 



27r2 
1 

1 



27r2 



3K 



h 



1 + K,^ 
3hv 



1 + 
^fX 



+ 



V>X 



9hl 



+ 



'Vx 



1 + 



+ 



Ah 



(75) 
(76) 
(77) 
(78) 
• (79) 



At this point one could easily determine the flxed points of the above system, 
investigate their stability, etc., but this is outside the scope of the present 
article. Rather, we note that since the x-field remains massive during the 
transitions studied in section ^ it decouples in the IR, in the sense that 
oo like l/A^ as A 0. For small A we are therefore left with the 
reduced system 



~~dt 

d\ 
"dt 



2., + ^^, 
^ 27r2 1 + ' 



k„ - 



9hl 



27r2 (1 + 



(80) 
(81) 
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This system has the same two fixed points we found in the Xcp^ case above, 
namely the trivial one and the WFP at k,^ = —1/7, = 87r^/49, = 
^vx ~ 0- What has changed is the stability of these points. While the 
Gaussian fixed point is still completely unstable, the WFP now has three 
unstable directions instead of just one. Thus it seems that fine-tuning one 
linear combination of UV couplings (by adjusting T) is not sufficient for 
a fiow to end up at the WFP. This, combined with the absence of other 
sufficiently stable fixed points, seems to indicate that the system does not 
undergo second order phase transitions. 

Given the above analysis one may well ask why the order parameter and 
the couplings seem to vanish continuously near the critical temperature, as 
shown in Figs. ^ and [^. The reason is that K^p and fiow to the WFP 
despite the instability in the other two directions in coupling space. This 
happens because the yj-sector of the theory completely decouples from the 
X-sector (but not vice versa), as can be seen from (p0|)-(p3|). In fact, (pO|) 
and (|8TD are exactly the same as the equations for the simple ip^ case, ([H]) 
and ([72|). Consequently all the reasoning for that case goes through and it 
suffices to tune one parameter (the temperature) in order to fiow to the fixed 
point|^. For our quartic truncation the transition is thus predicted to be 
second order and m^(Tc) and A^(Tc) fiow as in (|73 
behavior of h^p^ and at the phase transition can be found by plugging the 
fixed point values K^p = —1/7 and = 87r^/49 into ( ^2]) and (^3]). These 
equations are then easily integrated, with the result h^^ ~ exp(2t/3) and 

~ Ciexp(t) — C2exp(4t/3), where ci and C2 are positive constants. For 
A ^ we thus obtain X^^{A,Tc) ~ A^/^ and Ax(A,Tc) ~ ci - Ca/A^/^. This 

^•^It is actually not quite that simple. The fact that the WFP is a saddle point in the ip 
subspace indicates that one needs to fine-tune only one linear combination of couplings to 
flow to it. The important question is then for what values of the parameters adjusting the 
temperature enables one to achieve a correct linear combination. This is the question we 
have studied in detail in section ^, with the conclusion that the allowed parameter values 
are roughly those predicted by perturbation theory. In this appendix we argue only that 
if it is possible to produce a symmetry breaking phase transition by varying T, then this 
transition is second order and the critical behavior is equivalent to that of the Z2 model. 



and (0) as A ^ 0. The 



37 



behavior is observed in Fig. |^. The fact that A^^ does not go negative in the 
figure is simply due to the temperature resolution used. For T very close to 
one indeed finds A^(T) < 0. 
We conclude with a few comments. First, we have verified numerically 
that close to the RG trajectories of the full theory do get attracted to the 
WFP in the (p subspace. This shows that the WFP is the relevant one for our 
phase transitions. Second, we point out that the decoupling of the </3-sector 
from the x-sector at the phase transition will occur even if the evolution 
equations are not truncated at quartic order in the fields as was done above. 
The reason for this is simple: any diagram that contributes to the running 
of some coupling A^"^ (corresponding to an operator y?") is either built from 
vertices that have no x l^gs or else necessarily contains closed x-loops- The 
former are the same vertices present in theory, and the latter will be 
highly suppressed in the IR when = /i^/A^ oo. The critical behavior of 
our two scalar theory near the inverse phase transitions under consideration 
is thus equivalent to that of the Z2 model for any polynomial truncation 
of the effective potential. Finally, a comment regarding the behavior of A^ 
at Tc. While it may seem odd that this coupling goes to minus infinity as 
the cutoff is lowered, this does not mean that the theory is unbounded. It 
must be remembered that without our truncation higher order terms would 
be present and that the four dimensional couplings do not have a simple 
physical interpretation in the critical theory. To illustrate this point consider 
the Z2 model, and the coupling As corresponding to the operator ip^. At high 
temperature the dimensionless effective three dimensional coupling relevant 
for the investigation of fixed points is = T^AAg. The fixed point value 



of hs turns out to be negative [|T4|, which means that A8(A, Tc) —00 as 
A — *^ 0. In fact it is clear that all A„ with n > 8 will diverge at Tc. The 
point is that the physics near the phase transition is parameterized by the 
critical exponents of the effective three dimensional theory and not by the 
four dimensional couplings. 
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